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Abstract 

The adsorption of ethylene molecule on neutral, anionic and cationic gold clusters consist- 
ing of up to 10 atoms has been investigated using density-functional theory. It is demonstrated 
that C 2 H 4 can be adsorbed on small gold clusters in two different configurations, correspond- 
ing to the n- and di-a-bonded species. Adsorption in the % -bonded mode dominates over the 
di-a mode over all considered cluster sizes n, with the exception of the neutral C 2 H 4 — Au 5 sys- 
tem. A striking difference is found in the size-dependence of the adsorption energy of C 2 H 4 
bonded to the neutral gold clusters in the % and di-a configurations. The important role of 
the electronic shell effects in the di-a mode of ethylene adsorption on neutral gold clusters is 
demonstrated. It is shown that the interaction of C 2 H 4 with small gold clusters strongly de- 
pends on their charge. The typical shift in the vibrational frequencies of C 2 H 4 adsorbed in the 
%- and the di-a configurations gives a guidance to experimentally distinguish between the two 
modes of adsorption. 
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Introduction 



The adsorption of unsaturated hydrocarbons on transition metal surfaces has been studied exten- 
sively in order to understand the nature of hydrocarbon - metal interaction and chemical processes 
on solid surfaces; see, e.g., Refs.-^ 1 ^^ and references therein. The most significant attention 
was paid to investigation of ethylene adsorption,"2^»2J£JJ.iJ24MilliJ^Al^ilI because it is the simplest 
alkene containing an isolated carbon-carbon double bond. Hence it can be treated as a prototype to 
study the interaction and reactivity of different alkenes on metal surfaces. Moreover, the ethylene 
epoxidation is one of the most important processes in the chemical industry, because the product of 
such a reaction - ethylene oxide is widely used in various applications; see, e.g., Re.fs.±£»±M2*M12± 
and references therein. 

Ethylene can adsorb on metal surfaces in two different configurations .-LUiI 2 ^ 4 The first one is 
the it mode where one metal atom on the surface is involved in the adsorption of ethylene via a 
^-bonding. The second one is the di-c-bonded mode, when two metal atoms are involved in the 
adsorption via a o bonding. It was found that the di-a mode of adsorption is characterized by the 
increasingly important role of sp 3 hybridization in ethylene, while sp 2 hybridization (typical for 
free molecules) remains unchanged for ^-bonded species.— 

The bonding of ethylene with transition metals involves electron transfer from the filled bond- 
ing % orbital of the ethylene to the metal, alongside a back-donation from the d-orbital of transi- 
tion metal to the empty n* anti-bonding orbital of ethylene in accordance with the Dewar-Chatt- 
Duncanson model. 22 Thus the appearance of the sp 3 character of hybridization in the di-c-bonded 
ethylene can be explained by the increasing role of the electron back-donation from the metal to the 
Tt* anti-bonding orbital of ethylene. Hence the di-c-bonded ethylene is activated more strongly in 
comparison with the 7T -bonded one, thus it can be more reactive.-^ 2 - Therefore, an understanding 
of the specific mechanisms of ethylene adsorption on metal surfaces is important in order to gain a 
better insight of the catalytic processes and reactivity of adsorbed hydrocarbons. 

Surprisingly, not many works have been devoted to the investigation of ethylene adsorption on 
metal clusters and nanoparticles, despite the fact that chemical and physical properties of matter 
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at nanoscale are very different from those of the corresponding bulk solids. These properties are 
often controlled by quantum size effects and significantly depend on the size and structure of 
atomic clusters. 23 

A remarkable example is gold. It is well known that gold in its bulk form is a catalytically in- 
active and inert metal. However gold at nanoscale manifests extraordinary catalytic activity which 
increases with a decrease in the cluster size of up to 1-5 nm.— Moreover, recent studies 
demonstrate that catalytic activity of gold clusters adsorbed on an iron oxide support correlates 
with the presence of very small clusters of ~ 10 atoms.— It was also reported that gold clusters 
with the number of atoms n = 3 — 1 1 possess extraordinarily high electrocatalytic activity toward 
the 2 reduction reaction in acid solutions. 29 The origin of such size-dependent catalytic activity 
of gold remains highly debated and has yet to be fully understood. A comprehensive survey of the 
field can be found in review papers and books; see, e.g., R e f s 22a22i2L22i22aM it was demonstrated 
that the catalytic activity of gold clusters depends on the type of the support material, the pres- 
ence of defects in the material (e.g. F-center defects), the environment or special additives; see, 
e.g., R e fs,.2£i2L2£i2ai2Z and references therein. However, experiments indicate that gold clusters de- 
posited on an inert support can efficiently catalyze the process of styrene oxidation by dioxygen.— 
Hence, the origin of catalytic activity derives from the cluster itself and therefore even free clusters 
can be effective catalysts. 

Currently, a large variety of catalytic reactions on gold clusters has been studied experimen- 
tally. This includes the processes of catalytic oxidation of carbon monoxide, as well as more 
complex reactions such as alcohol oxidation, the direct synthesis of hydrogen peroxide and alkene 
epoxidation. 26 i 27 i 38 i 39 

However, theoretical studies of the gold nanocatalysis have mainly focused on the investigation 
of adsorption and catalytic reaction of 2 and CO.— ^ Nonetheless, the theoretical study of the 
interaction of alkenes with gold clusters has been relatively unexplored with the exception of works 
on the binding of propene on gold and mixed gold-silver clusters,— 1 ^ propene epoxidation 42 and 
our recent work on cooperative adsorption of oxygen and ethylene on small gold clusters. 43 



3 



In this paper we report the results of a systematic theoretical study of the adsorption of C 2 H 4 
on free neutral, anionic and cationic gold clusters consisting of up to 10 atoms. We find that C 2 H 4 
adsorbs on small gold clusters in two different configurations, corresponding to the %- and di-a- 
bonded species. Adsorption in the 7T -bonded mode dominates over di-a mode over all considered 
cluster sizes n with the exception of the neutral C 2 H 4 — Au 5 system, where the di-a-bonded con- 
figuration is energetically more favorable. We find a striking difference in the size dependence of 
the adsorption energy of C 2 H 4 bonded to neutral gold clusters in the % and di-a configurations. 
Thus, the adsorption energy of the ft -bonded C 2 H 4 develops non-monotonically as a function of 
cluster size with a local minima at n=6. The adsorption energy, calculated for the di-a-bonded 
C 2 H 4 exhibits pronounced odd-even oscillations, showing the importance of the electronic shell 
effects in the di-a mode of ethylene adsorption on gold clusters. 

We also demonstrate that the interaction of C 2 H 4 with small gold clusters strongly depends 
on their charge. The excess of a positive or a negative charge on the gold cluster can change the 
balance between electron donation and back-donation processes. Hence ethylene adsorption and 
reactivity can be manipulated by the cluster charge. This effect can be particularly important for 
understanding the mechanisms of catalytic activation of ethylene adsorbed on small gold clusters. 

Theoretical Methods 

The calculations are carried out using density-functional theory (DFT). The hybrid Becke-type 
three-parameter exchange density functional paired with the gradient-corrected Perdew-Wang 91 
correlation functional (B3PW91)— is used. The choice of the B3PW91 functional is consistent 
with our previous work on cooperative adsorption and catalytic oxidation of ethylene on small gold 
clusters. 43 The obtained theoretical data on dissociation energies and bonding in Au 2 (1.904 eV, 
2.546 A), and C 2 H 4 (7.596 eV, 1.324 A) are in a good agreement with those of earlier experimental 
studies, Au 2 (2.30 eV, 2.472 A),- 4 ^ and C 2 H 4 (7.76 eV, 1.339 A). 46 The standard LANL2DZ basis 
set of primitive Gaussians is used to expand the gold cluster orbitals formed by the 5s 2 5p 6 5d l0 6s l 
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outer electrons of Au (19 electrons per atom). The remaining 60 core electrons of the Au atom 
are represented by the Hay-Wadt effective core potential accounting for relativistic effects. 47 For 
carbon and hydrogen, the aug-cc-pVTZ basis set— is employed. Calculations have been carried 
out with the use of the Gaussian 03 code. 49 
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Figure 1: Optimized geometries of the most stable anionic, Au n , neutral, Au n , and cationic, Au„ , 
gold clusters calculated within the size range n= 1, ., 10 in the B3PW91/LanL2DZ approximation. 



The structural properties of neutral, anionic and cationic gold clusters have been the subject 
of numerous theoretical investigations; see, e.g., Refs.— >2L22 for a review. In the present work 
the cluster geometries have been determined with the use of the cluster fusion algorithm 50 which 
belongs to the class of genetic global optimization methods .■ 5 -^2 We have successfully used a 
similar approach to find the optimized geometries of various types of atomic clusters.— The 



optimized structures of the gold clusters Au n , Au n , and Au + (n = 1 — 10) obtained in the present 



work are presented in Figure 1 These structures are in a good agreement with those reported in 
previous theoretical studies; see, e.g., Refs.— i2k£L2!L32 

In order to obtain the most stable configuration of C 2 H 4 adsorbed on neutral, anionic and 
cationic gold clusters, we have created a large number of starting geometries by adding C 2 H 4 
molecule in different positions (up to 30) on the surface of the most stable cluster and up to eight 
isomer structures of the corresponding Au n , Au~, and Au* clusters. The starting structures have 
been optimized further without any geometry constraints. 

Theoretical results 

Geometry optimization for ethylene adsorbed on neutral, cationic and anionic 
gold clusters 

The results of the cluster geometry optimization for neutral C 2 H 4 — Au n , anionic C 2 H 4 — Au~ and 
cationic C 2 H 4 — Au + clusters within the size range n < 10 are shown in |Figure 2\ |Figure 3| and 
Figure 4[ respectively. 



Gold clusters with adsorbed ethylene possess various isomer forms whose number grows dra- 
matically with cluster size. Among the huge variety of different isomers, one can distinguish two 
specific adsorption modes that correspond to the %- and di-c configurations of the adsorbate. The 
appearance of two different modes of ethylene adsorption on small gold clusters is similar to that 
reported for ethylene adsorption on Au(100) and Au(l 1 1) surfaces. 14 However in the case of ad- 
sorption on the bulk, there are only two possible geometrical configurations of ethylene (we do 
not consider here \l -bridging species, 14 because ethylene adsorption in such a configuration is not 
possible in the case of small clusters). Moreover, there is no structural changes of the bulk surface 
due to adsorption of ethylene. In the case of ethylene adsorption on small gold clusters, each of 
the it- and di-c configurations possess a variety of structural isomers. These isomers correspond 
to the adsorption in nonequivalent positions on the surface of clusters. In |Figure^j |Figure3] and 
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Figure 4|we present the most stable structures of it- and di-a -bonded ethylene adsorbed on small 



gold clusters. The C - C and C - Au interatomic distances are given in angstroms. 

The interaction of gold clusters with ethylene results in an electron donation from the C 2 H 4 
to the gold cluster alongside a back-donation from the cluster to the adsorbate in accordance with 
the Dewar-Chatt-Duncanson model. 22 The rearrangement of the electronic structure and electron 
transfers upon adsorption lead to a dramatic change in cluster geometry. The interplay of electronic 
and geometry shell effects, also influence the cluster structure. 55 

Free neutral gold clusters possess planar, two-dimensional (2D) structures when the cluster size 
is up to 13 Au atoms. 59,60 In Figure 1[ we present the most stable structures obtained for Au n , Au + 



and Au" clusters (n = 1, 10) within the B3PW91/LANL2DZ DFT method. |Figure2] demon- 
strates that adsorption of C 2 H 4 in the n configuration leads to a prevalence of three-dimensional 
(3D) structures of gold clusters for n > 9; the structures of small Au 3 and Au 4 clusters transform 
to that which are energetically favorable for free cations Au 3 + and Au 4 + , respectively. For n < 9, 
ethylene favorably adsorbs at the cluster edge where the carbon-carbon bond lay in the plane of 
the cluster with the exception of the C 2 H 4 — Au 6 system, where the carbon-carbon bond is located 



perpendicular to the cluster plane (see Figure 2 ). 

Adsorption of C 2 H 4 in the di-a configuration results in a considerable rearrangement in the 
structure of neutral gold clusters. Thus, the geometry of Au 3 , and Au 6 - Au 10 clusters changes 
upon C 2 H 4 adsorption. It is known that gold clusters possess a large number of energetically 
close-lying structural isomers.— ^ Interaction with adsorbate can easily alter their energy order. 
It is important to note, that interaction with adsorbate can be responsible for structural rearrange- 
ments even for the relatively large clusters. Thus, recent DFT calculations demonstrate that the 
morphology of Au 79 nanoparticle can be transformed by interactions with a CO atmosphere. 61 

It is not surprising that the structures of positively and negatively charged ionic gold clusters 



change upon ethylene adsorption. Figure 3 demonstrates that although anionic gold clusters Au n ~ 
with n = 1, 10 remain planar upon ethylene adsorption (both for the % and di-a adsorption 
configurations), they manifest considerable structural rearrangements. Thus, gold cluster anions 
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a) 




Figure 2: Optimized geometries of neutral C 2 H 4 — Au n clusters: (a) 7T-mode of adsorption; (b) di- 
e-mode of adsorption (in the case of a single Au atom only one a bond is formed). The C - C and 
C - Au interatomic distances are given in angstroms. 
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Figure 3: Same as in Figure 2 but for anionic C 2 H 4 — Au n clusters. In the case of a Au 2 dimer 
only one o bond is formed. 
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Figure 4: Same as in|Figure"2| but for cationic C 2 H 4 — Au + clusters. 
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alter their structure upon ethylene adsorption in the ^-configuration for n = 4, 7, 8, and 10. In 
the case of di-a mode of ethylene adsorption, structural rearrangements occur for gold cluster 
anions with n = 3 - 8. A carbon-carbon bond of the adsorbed ethylene lie in a cluster plane for all 



structures presented in Figure 3 



Figure 4| shows that in the case of ethylene adsorption on gold cluster cations, Au+, the result- 
ing cluster systems C 2 H 4 — Au + exhibit 3D structures for n= 5 - 10 in the case of the % mode of 
ethylene adsorption; and n = 8 and 10 in the case of the di-a mode. In both cases of the n and di-a 
modes of ethylene adsorption, the structural rearrangements in gold cluster cations occur for n > 
5. 

Such a severe change in the morphology of gold clusters due to the interaction with ethylene 
makes the procedure of structural optimization rather complicated. In order to get energetically fa- 
vorable structures of C 2 H 4 — Au n , C 2 H 4 — Au n , and C 2 H 4 — Au* clusters, it is necessary to consider 
adsorption of ethylene not only on the most stable structures of the corresponding Au n , Au ~, and 
Au + clusters, but also on a large number of their energetically less favorable isomer states. In the 
present work we have taken into account up to 8 isomer structures of the host gold clusters. This 
insures that we do not miss energetically favorable structures with considerable change in cluster 
geometry due to C 2 H 4 adsorption. 

Energetics of ethylene adsorption on neutral, anionic and cationic gold clus- 
ters 

We first study the energetics of C 2 H 4 adsorption on neutral gold clusters Au„ with 1 < n < 10. 
The spin states of the optimized C 2 H 4 — Au n structures are doublet and singlet for odd and even n, 



respectively. Figure 5 shows the evolution of the molecular adsorption energy, £ads> calculated for 
the most bounded % and di-a configurations of adsorbed C 2 H 4 as a function of cluster size n. The 
adsorption energy is defined as 

£ ads = E m (Au„ ) + E tot (C 2 H 4 ) - E tot ( Au„ - C 2 H 4 ) , ( 1 ) 
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where E tot (M) denotes the total energy of the most stable structure of the molecule (cluster) "M". 




n 



Figure 5: Molecular adsorption energy, E & ^, calculated for C2H4 adsorbed in n (solid line) and 
di-cr (dashed line) configurations on the neutral Au n clusters with 1 < n < 10. Note, that in the 
case of a single Au atom only one a bond is formed. 

In our previous work, we have studied adsorption of ethylene on small neutral gold clusters 
in the 7T -bonded configuration. 43 Here, we demonstrate that C 2 H 4 can be adsorbed on small gold 



clusters in two different configurations, corresponding to the %- and di-a-bonded species. Figure 5 
shows that the 7T -bonded mode of C 2 H 4 adsorption is energetically favorable in the considered 
range of cluster sizes n, with the exception of a single Au atom which interacts with one carbon 
atom in C 2 H 4 , forming a a bond, and the Au 5 cluster, where the di-a-bonded configuration of 
C 2 H 4 — Au n is 0.05 eV more bounded compared with the 7T -bonded one. 



It is seen from Figure 5 that the adsorption energy of the ^-bonded C 2 H 4 on Au n clusters 
exhibit a maximum at n = 3 and 4 followed by a fast drop with a minimum at n = 6. For 6 < n < 10 
the adsorption energy of ^-bonded C 2 H 4 shows a weak oscillatory behavior as a function of n. 

The n-evolution of the adsorption energy calculated for the di-a-bonded C 2 H 4 is considerably 
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different from that of the 7T -bonded C 2 H 4 . |Figure 5 [ demonstrates that the adsorption energy of the 
di-a-bonded C 2 H 4 has an odd-even oscillatory behavior with a local maxima for odd numbers of 
Au atoms n = 3, 5, 7, and 9. For clusters with an even number of n, the adsorption energy E^ of 
the di-a-bonded C 2 H 4 is small. Hence, C 2 H 4 readily adsorbs in the di-a configuration — but only 
for neutral gold clusters with an odd number of Au atoms. 

It is interesting to note that in the case of C 2 H 4 adsorption on a Au(l 11) surface, the di-a mode 
( E ads° = 0-6502 eV) is energetically more favorable in comparison to the % mode (E^ ds = 0.1545 
eV).— A similar effect is observed for C 2 H 4 adsorbed on a Au(100) surface where theoretical 
values of E^\~° = 0.7790 eV and E£, = 0.2456 eV have been reported.— It is known that the type 
of C 2 H 4 bonding depends upon the surface structure. It has been found that the 7T -bonded mode 
dominates in the adsorption of low coordinated atoms and step sites. 62 In the case of small metal 
clusters, almost all atoms are essentially surface atoms that are undercoordinated as compared 
to the bulk. 63 Thus, the 7T -bonded configuration of the adsorbed C 2 H 4 dominates for small gold 
clusters. As the number of low coordinated atoms decreases with increase in cluster size, we 
suppose that the di-a mode will become dominant for large gold clusters. 



Figure 6 and Figure 7 present the n-evolution of the the C-C bond distance, rc-c> an d the total 



charge transfer to the adsorbate, Ag C2H4 , obtained from the natural bond orbitals (NBO) population 
analysis.— Figure 6 shows that the C-C bond in the ^-bonded C 2 H 4 is enlarged up to 1.37-1.39 
A in comparison with the case of a free C 2 H 4 (1.32 A). The tq-c calculated for the 7T -bonded 
C 2 H 4 depends slightly on the cluster size, evolving with n similar to £ads- On the other hand, for 
the di-a-bonded ethylene the C-C bond is enlarged up to 1.42-1.50 A, which is comparable with 
the length of a single C-C bond in ethane, C 2 H 6 . The n-dependence of tq-c f° r the di-a-bonded 
C 2 H 4 manifests an irregular behavior with maxima at n = 2, 4, 7 and 9. An increase in the C-C 
bond length, as well as as noticeable change of the bending of the H atoms for di-a-bonded C 2 H 4 
(see |Figure 2[ ), can be explained by the increasingly important role of the sp 3 hybridization due to 
the electron donation from the gold cluster to the anti -bonding n* orbital of C 2 H 4 . Therefore, rc-c 
in the di-a-bonded ethylene approaches its typical values for a single C-C bond in sp 3 hybridized 
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Figure 6: The C-C bond distance, rc-c, in C 2 H 4 adsorbed in the n (solid line) and di-a (dashed 
line) configurations on the neutral Au n clusters with 1 < n < 10. The dashed-dotted lines denote 
the equilibrium C-C bond distances for free C 2 H 4 and C 2 H 6 . Note, that in the case of a single Au 
atom only one a bond is formed. 

C 2 H 6 . Hence the di-a-bonded ethylene is activated strongly in comparison with the 7T -bonded one, 
thus it can be more reactive in catalytic processes. 

The total charge transfer to the adsorbed C 2 H 4 depends on the balance between the donation 
and back-donation processes. |Figure 7 demonstrates that in the case of the n mode of adsorption 
the total charge transfer Aq ClU4 is relatively small and evolves with n similar to yq-c- For the di-a 
mode of adsorption the back-donation from the gold to the n* anti-bonding orbital of ethylene 
dominates over electron transfer from C 2 H 4 to the cluster. The total charge transfer to C 2 H 4 
adsorbed in a di-a configuration is large and varies in the range of 0.29 - 0.43 a.u. for n = 2 - 9. 
It is seen from Figure 7 that Ag C2H4 calculated for the di-a-bonded C 2 H 4 exhibits a pronounced 
odd-even oscillations as a function of cluster size n with the exception of n =1. Note, that for n =1 
the charge transfer Aq Clli4 is small because in this case only one Au— C a bond is formed. Similar 
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Figure 7: The electron charge transfer, Ag C2H4 , to C 2 H 4 adsorbed in the it (solid line) and di-a 
(dashed line) configurations on the neutral Au n clusters with 1 < n < 10. Note, that in the case of 
a single Au atom only one o bond is formed. 

odd-even oscillations have been found in the n-dependence of the adsorption energy of 2 on Au n 
clusters; see, e.g.— i 5 - 6 -" 6 - 5 - and references therein. However, the adsorption energy of CO does not 
present any odd-even effects. 65 

The appearance of the odd-even oscillations in and Aq ClYl4 for the di-<7 mode of adsorption 
is a result of the electronic shell effects— in the gold clusters and can be described within the 
jellium model. In spite of its simplicity, the jellium model can explain, at least on the qualitative 
level, many physical properties and chemical reactivity of metal clusters.— Mi£2, Indeed, each Au 
atom has one 6s electron which delocalizes in the whole volume of the cluster. The delocalized 
electrons are moving in the field of a uniform positive charge background, which binds the valence 
electron cloud. The gold clusters with an even number of atoms (valence electrons) have closed 
electronic shell structure. Hence these clusters are more stable and less reactive in comparison 
with open shell clusters possessing an odd number of atoms (valence electrons). 
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The prevalence of the electron transfer from the cluster to the adsorbate (back-donation pro- 
cess) is responsible for the promotion of the di-a-mode of C 2 H 4 adsorption. Indeed, the di-a mode 
of C 2 H 4 adsorption on small gold clusters with an odd number of atoms leads to the transfer of 
unpaired electron from the gold cluster to the ethylene anti-bonding n* orbital, causing the weak- 
ening of the C-C bond. The neutral gold clusters with even number of valence electrons possess a 
closed electronic shell structure, hence electron transfer to the adsorbate and formation of the di-a 
bonded species is suppressed. 
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Figure 8: Molecular adsorption energy, £ads> (upper row); and C-C bond distances in adsorbed 
molecules, rc-c, (lower row); calculated for C 2 H 4 adsorbed on the anionic Au^ (left column) and 
cationic Au^ (right column) clusters with 1 < n < 10. Solid and dashed lines correspond to the % 
and di-a configurations of adsorbed C 2 H 4 , respectively. Dashed-dotted lines in lower row figures 
denote the equilibrium C-C bond distances for free C 2 H 4 and C 2 H 6 . 



Figure 8 demonstrates the evolution of molecular adsorption energy, E a ^ s and C-C bond dis- 



tances, rc-c, as a function of cluster size n for C 2 H 4 adsorbed on the anionic Au n (left column) 
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and cationic Au^ (right column) clusters with 1 < n < 10. 

The adsorption energy of C 2 H 4 on a gold cluster depends on the balance between donation 
and back-donations processes, hence it can be manipulated by the cluster charge. It is seen form 
Figure 8 1 and Figure 8b that excess of the negative or the positive charges on gold clusters strongly 



influence the molecular adsorption. Thus, the interaction of C 2 H 4 with the negatively charged gold 
clusters Au n ~ becomes weaker if compared with the neutrals; while excess of the positive charge 
results in strengthening of the C 2 H 4 and Au„ bond. 



Figure 8t i demonstrates that the ft-mode of C 2 H 4 adsorption is not stable for Au and Au 2 



The adsorption energy of the ^-bonded C 2 H 4 on the gold cluster anions Au~ is 0.17 eV in the 
case of a trimer Au 3 ~, and this increases with an increase in the cluster size: up to 0.82 eV at n =7 
and dropping to 0.54 eV at n =8, while slowly increasing with an increase in the cluster size up to 
n=10. 



Figure 8b shows that the size dependence of E ads for 7T -bonded C 2 H 4 on the gold cluster cations 



Au* is very different from that of Au ~ . Thus, C 2 H 4 readily adsorbs on the cationic monomer Au + 
and dimer Au 2 with an adsorption energy of 2.95 eV and 3.58 eV, respectively. Further increase 
in cluster size results in the sharp decrease in E a i s up to 2.06 eV for n =3 followed by its slow 
decrease for 3< n <6. At n =7 the adsorption energy E ads exhibits a maximum followed by the 
second sharp decrease at n = 8 and 9. 

For any n in the range 1< n <10, the adsorption energy of the 7T -bonded C 2 H 4 on cationic 
E£ ds (Au+), neutral, E£ ds (Au n ), and anionic, E^jAu~), gold clusters satisfies the relation: E^ ds (Au^) > 
EZ ds (Au n )>E* ds (Au-). 

Indeed, the excess of the positive charge on the cluster favours the process of electron transfer 
from the C 2 H 4 to the gold cluster, which stabilizes the it mode of adsorption, and strengthens the 
Au-C bond. A similar effect has been found for adsorption of CO on small anionic and cationic 
gold clusters.—^- It was found that the adsorption energy of CO is larger in the case of a cluster 
of cations and smaller in the case of anions. With the increasing cluster size, the differences in 
adsorption energies for different charge states become smaller.— 
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It is important to mention that strengthening the Au-C bond for C 2 H 4 adsorbed on the gold 
clusters with an excess of the positive charge is responsible for a cooperative effect in simultaneous 
adsorption of C 2 H 4 and 2 on neutral gold clusters. 43 Indeed, the adsorption of 2 on Au n results 
in an electron transfer from the gold cluster to the anti-bonding 2n* orbital of 2 . In this case, 
C 2 H 4 effectively adsorbs on the positively charged gold cluster. A similar effect has also been 
predicted for adsorption of propene on 2 — Au n clusters.— It was supposed that the binding of 
propene to 2 — Au n should be stronger than the binding to the Au n cluster due to the effective 
charge transfer from the gold cluster to the oxygen molecule. 40 

The excess of the positive or the negative charge on gold clusters results in a considerable 
change in the adsorption energy of the di-a-bonded C 2 H 4 . As has been discussed above, the 
adsorption energy of the di-a-bonded C 2 H 4 on small neutral gold clusters exhibit an odd-even 
oscillations as a function of the cluster size. These oscillations appear as a result of the electronic 
shell effects in gold clusters: Au n with odd n have one uncoupled electron that can easily transfer to 
the empty n* anti-bonding orbital of C 2 H 4 . Such an electron transfer favors the di-a configuration 
of the adsorbed C 2 H 4 . 



Figure 8 1 demonstrates that there is no odd-even oscillations in the adsorption energy of the 



di-a-bonded C 2 H 4 on the anionic Au~ clusters. In that case, the cluster always possesses a weakly 
bounded electron that it can transfer to the adsorbate. 

In the case of the cationic Au„ clusters, the electron back-donation process to C 2 H 4 is sup- 



pressed for any n. Figure 8 ~> shows that E^ s calculated for the di-a-bonded C 2 H 4 on Au„ exhibit 
two profound minima at n =3 and 9. In accordance with the spherical jellium model, Au^ and 
Aug clusters possess closed shell electronic structure with 2 and 8 valence electrons, respectively. 
Hence, Au^ and Au 9 + are chemically inactive, in particular in the processes where electron transfer 
from the cluster to the adsorbate is responsible for formation of the chemical bond. 



Figure 8 ; and Figure 8 1 demonstrate the n-evolution of the C-C bond distance, rc-c, m C 2 H 4 



adsorbed on Au n and Au„ clusters, respectively. The C-C bond in the 7T -bonded C 2 H 4 on the 
anionic and the cationic clusters is enlarged up to 1.36-1.40 A irrespective to cluster charge. This 
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is 0.04 - 0.08 A larger than that for the free C 2 H 4 (1.32 A). The di-a-bonded ethylene is activated 
strongly in the case of adsorption on the anionic Au n ~ clusters. Thus, in the di-a-bonded ethylene 
the C-C bond is enlarged up to 1.53-1.42 A for the anionic and only up to 1.48-1.37 A for the 
cationic gold clusters. Such an effect has a clear explanation. The excess of electrons on a cluster 
promotes back-donation of an electron from the cluster to the n* anti-bonding orbital of C 2 H 4 , 
hence promoting its activation. The excess of the positive charge on a cluster suppresses the back- 
donation process, thereby suppressing the C 2 H 4 activation. Thus, by changing the cluster charge 
one can manipulate the cluster reactivity. 

Vibrational spectra of ethylene adsorbed on small gold clusters 

Analysis of the vibrational spectrum of adsorbed C 2 H 4 directly reveals the mode of C 2 H 4 adsorp- 
tion and the state of hybridization. This is due to the sensitivity of the C-C stretching frequency 
to the hybridized state of carbon atoms.— In 1985 Stuve and Madix introduced the so-called no 
parameter for characterizing rehybridization in C 2 H 4 upon adsorption on metal surfaces. 11 The 
no parameter depends on the shift of the C-C stretching frequencies in C 2 H 4 and ranges from 
for the free C 2 H 4 , to 0.38 for K[(C 2 H 4 )PtCl 3 ] (Zeise's salt), and finally to 1 for C 2 H 4 Br 2 . Zeise's 
salt and C 2 H 4 Br 2 have been chosen by Stuve and Madix as models describing pure n- and di-c- 
bonding of C 2 H 4 , respectively. Thus the larger the shift in vibrational frequencies (compared to 
the free C 2 H 4 ), the greater the degree of rehybridization. 

A similar approach can be applied for C 2 H 4 adsorbed on metal clusters. ITable II presents the 
fundamental frequencies of free C 2 H 4 calculated in the harmonic approximation and obtained from 
experiment.— 

The obtained theoretical results are in a good agreement with experimental data. Some discrep- 
ancy, especially noticeable in the high frequency range, might result from anharmonic effects. The 
anharmonic corrections can be taken into account, for example by the method of n-mode coupling 
representations of the quadratic force field. 73 However such a consideration goes far beyond the 
scope of the present paper. In order to distinguish between n and di-c modes of C 2 H 4 adsorption 
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Table 1 : Fundamental frequencies of C 2 H 4 



Mode 


TT • 1 

Harmonic, cm 


t- • <- 72 -1 
Experiment, cm 


Vi {b 2u ) 


826 


826 


V 2 (b3u) 


975 


943 


V3 (big) 


981 


949 


V4 {a u ) 


1057 


1023 


V5 {b 3g ) 


1237 


1236 


v 6 (a g ) 


1 inn 


11/10 


v 7 (*i B ) 


1467 


1444 


v 8 (a g ) 


1693 


1623 


v 9 (&1«) 


3131 


2989 


Vio (a g ) 


3146 


3026 


Vn (b 3g ) 


3205 


3103 


Vl2 (&2«) 


3233 


3106 



one should to calculate the shift in the C-C stretching frequencies of the adsorbed C 2 H 4 with re- 
spect to the free molecule. In that case the anharmonic corrections to the vibrational energy levels 
can be partly canceled. Therefore in this paper, we perform analysis of the vibrational spectra of 
the adsorbed C 2 H 4 in the harmonic approximation. 

Following Ref., 11 we have selected two vibrational modes: and v$, corresponding to the 
direct C-C and the in-plane CH 2 scissoring modes of vibrations, respectively. Information on the 
shift of the vibrational frequencies of C 2 H 4 upon its adsorption can be obtained from infrared or 



Raman spectroscopy. Figure 9 presents the frequency dependence of the Raman scattering activity 



calculated for a C 2 H 4 molecule adsorbed in the % ( |Figure 9fr ) and di-a ( |Figure 9b ) configurations 



on the neutral Au n clusters with n = 1, 10. The stretching modes corresponding to the V(, 



and Vg vibrations in C 2 H 4 are marked by the vertical arrows. It is seen from Figure 9i that the 
spectral dependences of the Raman scattering activity calculated for the ^-bonded ethylene on 
Au n clusters are similar for all the cluster sizes n considered. The first maximum at 960 cm 1 
corresponds to the symmetric and asymmetric swinging of H atoms in the direction perpendicular 
to the plane of C 2 H 4 — similar to the v 2 and V3 vibrational modes in free C 2 H 4 . These vibrations 
remain unchanged in the ^-bonded C 2 H 4 . Two intensive maxima at 1280 cm -1 and 1560 cm -1 
correspond to the V(, and Vs vibrational modes in C 2 H 4 . These lines are red-shifted in comparison 
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Figure 9: Raman scattering activity of the neutral C 2 H 4 — Au n clusters: (a) ^-bonded C 2 H 4 ; (b) 
di-cr-bonded C 2 H 4 . Gaussian broadening of half-width 45 cm -1 has been used. The stretching 
modes corresponding to the V(, and Vs vibrations in C 2 H 4 are marked by the vertical arrows. Note, 
that in the case of a single Au atom, only one o bond is formed. 
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with the free C 2 H 4 . 
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Figure 10: Shift in the harmonic frequencies, and v$, calculated for the % (solid lines) and di-a 
(dashed lines) configurations of adsorbed C 2 H 4 on the neutral Au n clusters with 1 < n < 10. Note, 
that in the case of a single Au atom, only one a bond is formed. 



Solid lines in Figure 10 demonstrate shifts Av^ and Av^ in harmonic frequencies and Vs 
due to the it mode of C 2 H 4 adsorption on small neutral gold clusters. For the 7T -bonded C 2 H 4 , 
the scissoring mode Vs is shifted by 30-40 cm -1 towards low frequencies if compared with the 
V6 mode. Both the Av£ and Av^ shifts vary in the range of 58 - 142 cm -1 and exhibit similar 
dependence on the cluster size with local minima at n=6 and 8. The relatively small frequency 
shifts Av^ and AVg correspond to the values of the no parameter typical for the ^-bonded C 2 H 4 . 



Figure 9 demonstrates that the spectral behavior of the Raman scattering activity calculated for 
the di-a-bonded C 2 H 4 is different from that obtained for the ^-bonded C 2 H 4 . Thus, for cluster 
sizes n= 2, 4, 5, 7 and 9, a broad band appears in the frequency region of 900 - 1300 cm -1 . This 
band results from the strong mixing of the C-C stretching mode V(, with the hydrogen swinging (V2 
and V3) and twisting (V4) modes. The coupling of the V2, V3, V4, and Vf, vibrational modes occurs 
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as a result of considerable deformations in C 2 H 4 . These deformations include the bending of H 
atoms in respect to the plane of a free C 2 H 4 due to the sp 3 hybridization in di-<7 bonded C 2 H 4 , 
as well as the twisting of CH 2 groups in respect to each other along the C-C bond. Such twisting 
is especially noticeable for C 2 H 4 — Au 2 and C 2 H 4 — Au 4 clusters. The strong coupling between 
different vibrational modes makes identification of V(, and Vs vibrations to be rather difficult: the 
mixed modes form new levels, neither of which represent the true V(, and v$ states. 

Appearance of the sp 3 hybridization character in di-a bonded C 2 H 4 results in a considerable 
red-shift in the stretching frequencies. Such shifts AVg 1_CT and AVg 1_<T in harmonic frequencies 



and V8 are shown by the dashed lines in Figure 10 



It is seen that the Av^ 1-0 shift is considerably larger than the AVg 1_CT which can be explained 
by the strong coupling between the V2, V3, V4, and V(, vibrational levels. The Avf 1_tT and AVg 1_(T 
frequency shifts varies in the range of 200 - 370 cm -1 for n=2, ... ,10. These values correspond to 
the no parameter typical for the di-a-bonded C 2 H 4 . Thus, the analysis of the vibrational modes 
of C 2 H 4 adsorbed on gold clusters allows one to identify the specific modes of C 2 H 4 adsorption 
and the rate of rehybridization upon adsorption. 



Conclusions 

In summary, we have demonstrated that C 2 H 4 can be adsorbed on small gold clusters in two differ- 
ent configurations, corresponding to the it- and di-a-bonded species. Adsorption in the 7T -bonded 
mode dominates over the di-a mode over all considered cluster sizes n, with the exception of the 
neutral C 2 H 4 — Au 5 system. We found a striking difference in the size dependence of the adsorption 
energy of C 2 H 4 bonded to neutral gold clusters in the it and di-a configurations. The adsorption 
energy, calculated for the di-a-bonded C 2 H 4 exhibits pronounced odd-even oscillations, showing 
the importance of the electronic shell effects in the di-a mode of ethylene adsorption on gold clus- 
ters. The different rate of hybridization in the %- and di-a-bonded C 2 H 4 can be responsible for the 
different rate of catalytic activation and hence, different reactivity of adsorbed C 2 H 4 . 
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We have also demonstrated that the interaction of C 2 H 4 with small gold clusters strongly de- 
pends on the cluster charge. Hence ethylene adsorption and reactivity can be manipulated by the 
cluster charge. The strengthening of the ethylene - cluster interaction for gold clusters with an 
excess of the positive charge is responsible for the cooperative effect in simultaneous adsorption 
of C 2 H 4 and 2 on neutral gold clusters. This effect can play an important role in the mechanism 
of catalytic oxidation of alkenes by dioxygen on the surface of gold clusters. 

Finally we have shown that the analysis of the vibrational modes of C 2 H 4 adsorbed on gold 
clusters allows one to identify the specific modes of C 2 H 4 adsorption, the rate of rehybridization 
upon adsorption, and thus the C 2 H 4 reactivity. 

In the present work, we have considered adsorption of ethylene molecule onto small gold clus- 
ters consisting of up to 10 atoms. Many interesting problems beyond the scope of the present work 
arise when considering multiple molecular adsorption on gold nanoparticles of larger sizes up to 1- 
5 nm, where the strong dependence of the catalytic activity of gold nanoparticles has been observed 
experimentally. In particular, multiple molecular adsorption on the surface of gold nanoparticles 
can result in a considerable change in cluster morphology, even for relatively large cluster sizes. 

Another important direction for future development is to investigate the effect of alloying on the 
adsorption and catalytic activation of alkenes on gold clusters. The doping of a gold cluster with 
atomic impurity can result in considerable changes in the geometrical and electronic structure of 
the cluster, thereby modifying and controlling the catalytic activity. Understanding how to enhance 
and control chemical reactions on a cluster surface is therefore a vital task for nanocatalysis. 
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